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In higher Landau levels (N > 0) and around filling factors v = 4/V + l, a two-dimensional electron 
gas in a double-quantum-well system supports a stripe groundstate in which the electron density in 
each well is spatially modulated. When a parallel magnetic field is added in the plane of the wells, 
tunneling between the wells acts as a spatially rotating effective Zeeman field coupled to the "pseu- 
dospins" describing the well index of the electron states. For small parallel fields, these pseudospins 
follow this rotation, but at larger fields they do not, and a commensurate-incommensurate transi- 
tion results. Working in the Hartree-Fock approximation, we show that the combination of stripes 
and commensuration in this system leads to a very rich phase diagram. The parallel magnetic field 
is responsible for oscillations in the tunneling matrix element that induce a complex sequence of 
transitions between commensurate and incommensurate liquid or stripe states. The homogeneous 
' and stripe states we find can be distinguished by their collective excitations and tunneling I — V, 

^ , which we compute within the time-dependent Hartree-Fock approximation. 
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PACS numbers: 73.43.-f, 73.21.-b 



I. INTRODUCTION 



'-^J \ In strong magnetic fields, the two-dimensional electron gas (2DEG) in a double quantum well system (DQWS) is 
£5 ■ well known to support spontaneous interlayer coherence if the interwell separation is smaller than some critical value 
d c . At filling factor v = 1, the interlayer .coherence is reponsible for many interesting properties of the 2DEG such as 
i a quantum Hall state, a Goldstone modea (in the absence of tunneling) associated with the broken U(l) symmetry of 

the pseudospin order in the ground state, which has been recently detected experimentallyEl, and exotic topological 
quasiparticle excitations such as merons and bimerons. When the magnetic field is inclined towards the plane of the 
2DEG in the v = 1 coherent state, a radical change in the behavior of the transport gap was observedB and explained 
as a commensurate-incommensurate transition in the pseudospin ordenj. A review of these properties can be found 
in Ref. | 



At bigger filling factors, in Landau levels N > 0, it has been shown! that the 2DEG in a DQWS can support states 
where both the electronic densities and the interlayer coherence are spatially modulated. Ref. || discussed the ground 
£N) ' state of the DQWS at filling factors v = AN + 1 when both wells have exactly the same number of electrons using the 
Hartree-Fock approximation. They showed that, for some range of interwell separations, a "staggered" stripe state 
is lowest in energy. This state involves a unidirectional charge density wave ordering of electrons in each well, such 
that the regions of high electron density in one well are nearest to regions of low electron density in the other. Stripe 
states in single well systems were first discussed by Koulakov, Fogler, and Shklovskii and Moessner and ChalkerEI and 
are believed to explain a large anisotropy observed in transport experiments for quantum HalLsystems in moderate 
magnetic fields, for which the highest energy occupied Landau level is approximately half-filledo. 

In a DQWS, stripe states are more complex and interesting than in the single layer case because of the possibility 
that interlayer coherence may coexist with charge ordering. In the staggered stripe state, interlayer coherence is 
maintained only along linear regions at the edges of the stripes. These "linear coherent regions" (LCR's), whose 
width decreases .with increasing interlayer separation, support topological quasiparticle excitations closely akin to sine- 
Gordon solitonsQ. In the absence of tunneling, Cote and Fertigy showed that this staggered stripe state supports two 
Goldstone modes associated respectively with broken translational symmetry and broken U(l) symmetry associated 
. - - i with interlayer coherence. 

In this paper, we address the questions of how the staggered stripe state and its collective excitations are modified 
by the presence of a magnetic field [Component parallel to the plane of the layers and how these changes are reflected 
in the tunneling current. We recalltj that, in the presence of a parallel magnetic field, the tunneling hamiltonian can 
be written as an effective Zeeman coupling of the form 



H 7 



fdr h(r)-S(r), (1) 



where S (r) is the pseudospin density and h (r) is a pseudo-magnetic field that rotates in the xy plane. In this paper, 
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we choose By along the y axis and use a gauge in which this pseudo-field is given by 

h (r) = 2t N (Q) [cos (Qx) x - sin (Qx) y] , (2) 

where Q = d/£p with d the separation between the wells, and £\\ = ^JhcjeB^ the magnetic length due to the in-plane 
component Sii of the magnetic field. At filling factor v = 1, the ground state at small layer separation and B\\ = 
is well approximated by a state in which all pseudospins are oriented along the x axis (we assume that the bare 
tunneling energy to is not zero). When the parallel magnetic field is turned on, the pseudospins rotate in the xy 
plane, aligning with the local orientation of h(r). This defines a uniform commensurate state, which we will call the 
commensurate liquid (CL). When Bu increases, the rotation period of the pseudospins decreases. This rotation is 
opposed by the Coulomb exchange interaction that tends to keep peudospins parallel. Above some critical Bu , there 
is a phase transition to an incommensurate state where the pseudospins no longer rotate with the in-plane field but 
instead tend to align along some common direction as Bu increases above the critical value. We call this state the 
incommensurate liquid (IL) state. It is this change in the ground state of the 2DEG with B\\ that is believed to be 
responsible for the radical change in the behavior of the transport gap observed by Murphy et al.u 

In the experiment of Murphy et at, the system consisted of two homogeneous two-dimensional electron gas (2DEG) 
with filling factor (in each well) v = 1/2 and in Landau level N = 0. The same commensurate- incommensurate 
transition should occur in higher Landau levels at total filling factors v = AN +1. In this case, however, we must 
recognize that, for d > d c (N), the system is unstable to the formation of stripes as was shown in Ref. In the 
pseudospin language, this staggered stripe state consists of pseudospin modulation in the xz plane. Applying a 
parallel magnetic field will tend to induce a commensurate rotation of the pseudospins in the xy plane in addition to 
the xz modulations. In this paper, we call this phase a commensurate stripe phase (CSP). If the parallel field is too 
strong, the exchange energy cost of this commensurate rotation will be prohibitive and the ground state will revert 
to a state with xz modulations only as if there were no tunneling and no parallel field. We call this other state an 
incommensurate stripe phase (ISP). As d and B\\ are changed for a given tunneling parameter, we thus expect a rich 
phase diagram with transitions between the four possible states just introduced: CL, IL, CSP and ISP. (We cannot, 
of course, rule out the possibility of other ground states with lower energy). 

In higher Landau levels, the effect of B\\ is more complex than for N = 0. This is because the tunneling amplitude 

in the iVtri Landau level has a non trivial B\\ dependence, t/v (Q) = toe^^^^L^ 2 ^ ) > w ith to the bare tunneling 
amplitude. For N > 0, this tunneling amplitude oscillates with wavevector Q = d/t», becoming negative in some 
range of B\\ . One effect of these oscillations is that B\\ can induce not only one but several (reentrant) commensurate- 
incommensurate transitions in the liquid or in the stripe phases depending on the value of N. Figure 1 illustrates some 
typical phase diagrams we find for several values of the bare tunneling parameter t$. These figures will be explained 
in more detail below. 

Another effect of the parallel field is that, in the commensurate stripe phase (CSP), the long-wavelength dispersion 
relations of the low-energy phonon and pseudospin wave modes will be strongly dependent on B\\ through ijy (Q)- 
We find that the gap in the pseudospin wave mode (which is related to the broken U(l) symmetry of the pseudospin 
order) and in the phonon mode (at finite wavevector) follows the oscillations of the tunneling amplitude t/v (Q)- In 
particular, the phonon and pseudospin wave modes become essentially gapless whenever the parallel field is such 
that t jy (Q) is close to one of its N zeroes. Such behavior is not seen in the incommensurate stripe phase (ISP) 
since the properties of this state are independent of t^ (Q) and so of Sm . We thus find that the CSP and ISP are 
clearly distinguishable by their low-energy collective modes. Figure 2 illustrates these modes as computed using the 
time-dependent Hartree-Fock approximation, as explained in more detail below. 

OnCjHiay to probe_these states is via interlayer tunneling. For very small values of to, the tunneling I — V has been 
arguedEj and showro to probe the Goldstone mode dispersion of the (spontaneously) coherent state as the parallel 
field is adjusted. An analogous effect occurs in the stripe state of this system, as illustrated in Fig. 3. A clear signature 
of the existence of stripes can be found by noting the non-monotonic behavior of the tunneling peak as a function 
of B11. In particular, one finds this peak at zero bias voltage for Bu = nhc/ade where a is the separation between 

the stripes in a given well and n = 0, 1,2, This is a reflection of the periodicity of the collective mode spectrum 

illustrated in Fig. 2(a), and directly demonstrates the existence of stripes in the groundstate. 

This paper is organized as follow. In Sec. II, we derive the tunneling Hamiltonian in the presence of the parallel 
magnetic field and explain the origin of the oscillations in t/v (Q). In Sec. Ill, we derive the Hartree-Fock equations 
of motion for the single-particle Green's function from which we compute the density and pseudospin density as well 
as the energy of the commensurate and incommensurate stripe states. The energy of these phases are compared with 
those of the commensurate and incommensurate liquid phases in Sec. IV to obtain the phase diagram, in Landau 
level N = 2, for various parameters Bu, to, and d. We briefly review, in Sec. V, the formalism necessary to compute 
the susceptibilities and dispersion relations in the time-dependent Hartree-Fock approximation (TDHFA) and present 
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FIG. 1: Phase diagram of the 2DEG in Landau level N = 2 in the B\\ — d/£± plane for values of the bare tunneling parameters 
equal to (a) to/ (e 2 / ' nt) = 0.01, (b) to/ (e 2 / nt) = 0.05 and (c) to/ (e 2 / 1 k£) = 0.1. The areas in white(gray) represent the 
commensurate(incommensurate) liquid state. The heavy line marked d c indicates the critical spacing above which one of the 
stripe state (either commensurate or incommensurate) is lower in energy than both the CL and IL phases. 





our numerical results in various region of the phase diagram studied in Sec. IV. Section VI shows how the tunneling 
current, in the stripe states, is affected by the parallel magnetic field. We conclude in Sec. VII. 



II. TUNNELING HAMILTONIAN 



We consider an unbiased symmetric DQWS in a tilted magnetic field B = £>||y + B±z at total filling factor 
v = 4iV+ 1 where N = 0, 1, 2, ... is the Landau level index. To simplify the analysis, we make several approximations. 
We retain only one electric subband of the DQWS, and assume that g*/i B B » A$as where Asas is the symmetric 
to antisymmetric gap so that the real spins are completely frozen. Furthermore, we consider the lower Landau levels 
as filled and inert and neglect any Landau level mixing. Considering very narrow wells and treating the tunneling 
term in a tight-binding approximation, the non-interacting Hamiltonian of the 2DEG can be written, in the gauge 
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FIG. 3: Variation of the tunneling current / (£) with the potential biais eV for several values of the parallel magnetic field 
in the CSP. The tunneling current is calculated by using the first term in Eq. @ and with TV = 2, to/ (e 2 / 'nl) = O.Of and 
d/£± — 1 .4. Note that Bu = 0.325 T corresponds approximately to Q£± = 27r/a where a is the separation between the stripes 
in a given well. 



where A = (0, B±x, — B\\x) , as 

H Q = Y.EoclcjCxj -^trt (Q) (e^ x c\ R c x ,L + e~ lQX c\ l c x ,r 



(3) 



where Eq is the electron energy in each isolated well and x is the creation operator for an electron in Landau level 
N, guiding center index X , and well index j = R, L. The electronic wavefunctions are given by 



N,X,j 



( r ) 1 e -«^ 2 ^ (x-X) Xi (z), 



(4) 



where Xj( z ) 1S the envelope wave function of the lowest-energy electric subband centered on the right or left well 
and ip N (x) is an eigenfunction of the one-dimensional harmonic oscillator. In the tight-binding approximation, the 
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FIG. 4: Variation of the tunneling amplitude tjv (Q) with the parallel magnetic field for filling fractions v = 1, 5, 9 corresponding 
to Landau levels TV = 0, f , 2 and for d/£± = f .0. 

tunneling amplitude, tpj (Q), is given by 

/+oo /(1 2 £ 2 \ 

dx^ N (x - X) e iQx <p N (x-X) = t e- Q e ^ 4 L% ( , (5) 

where to is the bare tunneling amplitude (defined as a positive quantity), L° N (x) is a generalized Laguerre's polynomial 
and the wavevector and Q = d/£ 2 where £\\ = ^Jhc/eB^ is the magnetic length associated with the parallel magnetic 
field. 

If N >-p, the tunneling amplitude oscillates with the strength of the parallel field and/or the separation between 
the wellalJ. (In this paper, we neglect the further dependency of to on d due to the change in the overlap of the electric 
subband wavefunctions Xj( z ) °f the two wells). We remark that the first zero of ijv (Q) appears at relatively small 
parallel field. For example, since 2irno£ 2 L = v where Uq is the total areal density of the electrons in the DQWS, the 
wavevector Q£± = v (d/£±) (e/n hc) Bn, and so d/£± = dy / 2irno/i/. If we assume a typical total electronic density 
of no = 1 x 10 11 cm~ 2 , we have, with Bii in Tesla, Q£± — 0.242 (y) {d/l±)B]\, On the other hand, at this density 
the separation between the wells can be expressed as (d/£x) = 0.793d(100A)/ v / i/. For total filling factor v = 1, we 
get d/£± = 1.6 for a typical interwell separation of 200 A. This value decreases as N increases at fixed total density. 
Figure 4 shows the behavior of ijv (Q) for d/£± = 1.0, n = 1 x 10 11 cm -2 and filling fractions v = 1,5,9 corresponding 
to N = 0, 1, 2. We see from this figure that for N = 2, t^(Q) vanishes near B\\ ss 0.5 T and 1.2 T, and that a parallel 
field of sa 3 T is required to reach the asymptotic behavior of the tunneling amplitude. 



III. HARTREE-FOCK EQUATIONS OF MOTION 

We describe the various ordered coherent states of the electrons in the partially filled Landau level N by the set of 
average values where ftj(q) is an operator that we defindlj by 

ftjW = ^ E e-^H**^ cl xCj , x _ qy e, (6) 
x 

where N v is the Landau level degeneracy. At a semiclassical level, the (/^(q)) s can be thought of as the Fourier 
transform of the density of guiding centers of the cyclotron orbits while the off diagonal terms (p^j(q)) describe 
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quantum coherence between the two 2DEG. The average values (p, j (q)) are obtained by computing the single-particle 
Green's function 

G itj (X, X\ r) = - (Tc liX (t) c\ x , (0)) , (7) 

whose Fourier transform we define as 

Gi,jM = ^Yl e- L M x + x ')5 x , x ,_ qvl *G i , j {X,X l ,T), (8) 
* x,x> 

so that Gij (q,r = 0") = (p^ (q)). 

Using the tunneling Hamiltonian Ho [Eq. (^)] and working in the Hartree-Fock approximation, we find that the 
single-particle Green's functions for the DQWS obey the system of equations 

[ituj n - (-Bo - a 1 )] Gr,r (q, w «) + ** (q) Gl,.r (q - Q, w «) (9) 

- ^2 Ur < r ( q > q/ ) Gr ' r ( q/j _ X! c/r » x ( q ' q ') ( q/ ' w ™) = fi(5 i^ ' 

q' q' 



[tfiw„ - (£ - aO] G fl , L (q, w„) + t* (q) G L , L (q - Q, u n ) (10) 

( q> q/ ) Gr > l ( q/ ' w ») ~ X! ( q ' q/ ) Gl > l ( q/ > = °' 



q' 



[ifiu n - (Eq - (i)] G l ,_r (q, uj n ) + t (q) G;^ (q + Q, u) n ) (11) 

- J! U l,R (q, q') Gfl, fl (q', c^ n ) - ^ L/^l (q, q') G LM (q', w n ) = 0, 
q' q' 

[ifiw„ - (E - /i)] G l ,l (q, w n ) + * (q) Gfl^ (q + Q, w n ) (12) 

- X! ^.-r (q> q') Grl (q'' w ») _ X! ( q ' q/ ) Gl ' l ( q '' = fi(5 i>°- 

q' q' 
The tunneling term t (q) is, in our particular geometry where the parallel magnetic field is along the y axis, given by 

*(q) =MQ)e 4W ^ 2 . (13) 

When Q = 0, Eqs. (p|)-(|l2]) is the system of equations we studied in Ref. || The Coulomb interactions and averaged 
densities in this system of equations produce the mean field potentials 



Ur,r (q q') = ( -j- J [H (q - q') - X (q - q')] ( Pli . R (q - q')) 7 q , q ' + ( ) H (q - q') (p ij£ (q - q')> 7 q , q 



C/l,l (q, q') = ( [H (q- q') - X (q - q')] (q - q')> 7 q , q ' + f ^ ) 5 (q - q') <p fl , fl (q - q')> 7,,,' 



(q q') = - X (q - q') (p L>R (q - q')) 7 q , q <, (14) 

U'i.fl (q q') = - (^-) X (q - q') (p flii (q - q')) 7 q , q ', (15) 



where 

7q = e -<qxq')' 2; i/ 2 . (16) 
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The definitions of the coulomb interactions H, X, H and X were given in Ref. Q They are respectively the Fourier 
transforms of the intrawell Hartree and Fock interactions and the interwell Hartree and Fock interactions. 

Once the densities are extracted from the single-particle Green's functions, the Hartree-Fock energy per electron 
can be written as 



E 



HF 



-t N (Q)Re[(p R r (Q)>] 



(17) 





( e 2 N 




(pfl,fi(q)>| ^ 


-\(pl,l(<D)\ 2 ' 




f e 2 N 


)£*(q) 

a 


\(Pr,r(<1)) \ 2 ^ 


-\(pl,l (q)>| 2 " 



e 

kI_i 
kI 



£#(q) [(pr,r{-^)(pl,l (q))] 

q#0 



^X(q)|( fe)fi (q)>| 



The (p 4j - (q))'s are found by solving Eqs. (10)-(13), and using Eqs. (7)-(9) to extract them from the resulting 

Green's functions. Since the potentials entering the former set of equations involve the (Pi j (q))'s, this procedure 
must be iterated until convergence is achieved. (The numerical scheme is described in more detail in Ref. E2I) We 
next describe the results of such calculations for several possible states, and discuss the phase diagram of the system. 



IV. NON-HOMOGENEOUS HARTREE-FOCK STATES IN HIGHER LANDAU LEVELS 

In Refs. H and [|, it was found, for Bu = 0, that, with the exception of the N = Landau level, the 2DEG ground 
state first evolves from a uniform coherent state (UCS) to a unidirectional CDW state, and then undergoes a first 
order transition to a modulated stripe phase as d/£± is increased from zero. Wigner crystal (WC) states were found 
to be higher in energy than the unidirectional modulated states for N > 0. In the N — Landau level, the evolution 
was shown to be from a uniform coherent state to a CDW state and then to a Wigner crystal state at large interwell 
separation. The coherent states (either CDW or WC) are always lower in energy than their incoherent counterparts. 

In this paper, we consider the consequences to the stability of these unidirectionally ordered states of adding a 
parallel magnetic field. In view of our earlier results, we restrict ourselves to higher Landau levels N > and consider 
only coherent states. (Incoherent states are insensitive to the parallel field). Our filling factor is v = 4N + 1 and we 
consider four possible groundstates. These states are best described in a pseudospin language where the right (left) 
well electronic state is associated with a pseudospin up (down) state. The p i j operators defined above can be mapped 
into the density and pseudospin density operators by the relations 



n(q) 


= 2 \. p R,R (q ) 


+- pl,l (q)] , 


(18) 


p x (q) 




-i- plm. (q)] , 


(19) 


p y (q) 


= \ Wl (q) 


- pl,r (q)] , 


(20) 


Pz (q) 


= 2 [ p «>-R ( q ) 


- pl,l (q)] • 


(21) 



The tunneling hamiltonian of Eq. (|3|) can be written as an effective "Zeeman" coupling 

H T = -Y,h(X)-p(X), (22) 
x 

where p (X) and n (X) are the (guiding center-) spin density and density operators respectively. The field h (X) plays 
the role of a pseudo-magnetic field that lies in the plane of the 2DEG and rotates in the space of the X coordinate 
according to 



h (X) = 2t N (Q) [cos (QX) x - sin (QX) y] , 



(23) 
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|23|) is due to our particular convention for associating 



with a period A = 2-k/Q = 2irlyd. (The minus sign in Eq. 
the R and L wells with the up and down spins). 

When the magnetic field is tilted towards the plane of the 2DEG (at constant B± to keep the filling fraction 
constant), the rotation period A of the pseudospins decreases. In a non-interacting picture, the pseudospins would 
align with the local field h (X) to minimize the tunneling energy. When electron-electron, interactions are taken 
into account, however, the pseudospins interact through an effective exchange intcractiona that favors a parallel 
configuration of the pseudospins. At small parallel magnetic field, the pseudospins can follow the effective Zeeman 
field. The resulting state is called the commensurate liquid (CL). Above some critical parallel magnetic field, Bn c , the 
exchange energy cost of the rotation is too high, and an incommensurate liquid (IL) state where the pseudospins are 
spatially uniform becomes lower in energy. We note that this incommensurate state is never the true ground state: at 
a field slightly lower than Bi\ c , the 2DEG has a transition into a soliton lattice state (SLS) where kinks (slips of 2ir 
in the phase of the rotating pseudospins) form a lattice in the direction perpendicular to the parallel magnetic field 
with a period X s (Q). This SLS continuously evolves into the IL as the parallel field is increased. We do not explicitly 
treat the SLS in our work since it is significantly different than the IL only in a small range of Q (see, for example, 
Rcf . |l3|) . We will thus refer to the commensurate-incommensurate transition as to the transition between the CL and 
the IL. 

The CL and IL states are spatially homogeneous; their orderings are characterized by the behavior of (/^^(q)) 
but have no non- vanishing values of {pn{(\ 7^ 0)). In high Landau levels, stripe states may fornoE] which have guiding 
center densities oscillating spatially between the wells. It was found in Ref. ^ that for d/£± not too small, the edges of 
the stripes in each well will admix to form linear coherent regions (LCR's), and that in the pseudospin language such 
states are akin to helical ferromagnetsa. In the presence of a parallel field, the pseudospin ordering in the LCR's may 
or may not follow the effective Zeeman field, leading to a commensurate stripe phase (CSP) and an incommensurate 
stripe phase (ISP) . 



A. Commensurate Liquid (CL) 



The commensurate liquid at v = 1 and with tjq (Q) > is defined by homogeneous and equal densities of electrons 
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in both wells and by a rotation of the pseudospins in the xy plane with a single period A = 2ir£f,/d. The order 



parameters of this state are given by 



(pr,r(°)) = (pl,l(P)) = -, 
(pr,l(-Q)) = (pl,r(Q)) = 1- 



(24) 
(25) 



All other p's vanish in this state. The energy per particle Ecl = Ehf/N c (N e is the number of electrons in the 
partially filled level) is simply 



EcL=-\t N {Q)\-\(^ r - ) V(()) + ViQ) 



(26) 



The absolute value in Eq. (g6|) is necessary since, when (Q) < 0, the pseudo-magnetic field reverses direction and 
we must have (p R L (— Q)) = — i. 



B. Incommensurate Liquid (IC) 

The incommensurate liquid is a state with equal and homogeneous densities in both wells and with the pseudospins 
all aligned in the same direction which we take here as the x axis. The only non-vanishing order parameters are 



(pr,l (0)),(pl,r(°)) 
(Pr,r(0)),(Pl,l (0)) 



and the Hartree-Fock energy per particle is 



Eic — — 



I 



X (0) + X (0) 



(27) 
(28) 

(29) 
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which does not depend on the tunneling term or the parallel magnetic field. 



C. Commensurate Stripe Phase (CSP) 

In what follows, the parallel magnetic field is assumed to be in the y direction so that Q = QqSc. The coherent 
stripes are defined by charge and spin modulations with reciprocal lattice vectors G = — ne where n = 0,±1,±2, ... 
and e is a unit vector in a direction making an angle a with Q and the x axis. In the absence of the parallel magnetic 
field, the coherent stripes would be described by the order parameters {\Pij (G))}. With the field, the situation is 
slightly more complex since the tunneling term forces (p R L (Q)) 7^ and in general Q ^ {G}. 

To apply our equation of motion method, we need a single set of wavevectors to describe the modulated state. This 
can be achieved by considering a solution for the transverse part of the spins, of the form 

(PR,L(r)) = e-^(p R . L (r)) : (30) 

where (p R L (r)\ has the periodicity of the stripe state [i.e., (p R L (r)\ 7^ for q G {G}]. The sole effect of the rotating 
field is to force a rotation of the pseudospins in the xy plane in the stripe state which itself has the pseudospins rotating 
in the xz plane. We call this solution the commensurate stripe phase (CSP). Finding the order parameters is most 
easily accomplished by going into a spatially rotating frame where the xy pseudospins are aligned. We define the 
average values of the operators in the rotating frame by 

(pr,l(<i)) = <pfl,i(q-Q)>, (3i) 
<Px,fl(q)) = <p^(q+Q)>- (32) 

Because the system of equations (||)-(|l^) contains the factor 7 q q / , making the substitution q — > q ± Q introduces 
additional phase factors in the equations. We can preserve the structure of this system of equations in the rotating 
frame if we also redefine the diagonal order parameters as 

(p£ fl (q)) = 7 *(q)(p^(q)), (33) 

(p^(q)) = 7 (q)(^(q)>, (34) 

where 

7 (q) = e iqxQi l /2 = e -*«v«£/2. (35) 



The CSP thus occurs when (q)y are non- vanishing for q e {G}. 

To carry out the computation, it is convenient to define a 2 x 2 Green's function matrix, with matrix elements 

G R , L (q,uj n ) = G R . L (q+ Q, w n ) , (36) 

Gi,B(q,w n ) ee G h ,R (q - Q, w„) , (37) 

GR, R (q,u; n ) ee 7* (q) G R . R (q,w„) , (38) 

G L:L (q,uj n ) ee 7 (q) G LtL (q, u) n ) . (39) 



Using the fact that Q = Qx, we have 



i(q±Q)=t(q). (40) 



With these definitions and identities, the system of equations for the single-particle Green's function in the rotating 
frame is 



1 , 

w„ - t ( E o - N 



Gr,r (q, u„) + i« (Q) Gl,r (q, w n ) (41) 



^ Ur,r (q, q') Gr.r (q', u) n ) - Ur,l (q, q') G l ,r (q', w n ) = £q,o, 

q' q' 
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iuj n - - (E - fj) 



Gr,l (q, w B ) + ~zt (Q) G L . L (q, UJ n ) 



q' q' 



iu n - 4 (-Bo - M) 



G l ,_r (q, w„) + h; (Q) Gr,r (q, w„) 



E ^.-R ( q ' q ') Gr.r (q', w n ) - C/l : l (q, q') G l ,_r, (q', w„) = 0, 
q' q' 



(42) 



(43) 



iuj n - - (E - (i) 
n 



G l ,l (q, w n ) + ~t (Q) (q, uj t , 
n 



- E Ul > r ( q > G ^ ( q '' w «) ~ X! ( q > q/ ) Gl > l ( q '' w «) = s ^°- 

q' q' 
The four mean-field potentials, in this frame, are given by 

Urm (q, q') = 7* (q) Ur,r (q, q') 7 (q') , 
Ul.l (q,q')= 7 (q) ^ (q, q') 7* (q') , 



(44) 



(45) 
(46) 



%,L(q,q') = %,L(q + Q,q')7*(q'), 
Ul,rM) = U LtR (q-Q,q') 7 (q') . 



Finally, the Hartree-Fock energy per particle is given by 

2 



CSP 



J2t N (Q)Be[(p R>L (0))] 



1 



/ e 2 N 


)E^( q ) 


|<^,i?.(q)> 2 + 


-\(pl,l (q))| 2 


q^O 






/ e 2 N 


)E^( q ) 

n 


IWi?.(q)> 2 -+ 





+- 



e 



E^( q ) i-f (q)] 2 [(^(-q))(PL,L (q)>] 



q#0 



K^)E^( q + Q)IK,(q)>| 



(47) 
(48) 



(49) 



At fixed d/l± and Bu, the energy of the CSP is obtained by minimizing the Hartree-Fock energy of Eq. (^) with 
respect to a, the separation between the stripes in a given well (which may be written in unitless form as £ = a/l±) 
and with respect to the orientation a of the stripes with the x axis. When a^O, the commensurate rotation of the 
pseudospins will induce a modulation in the total density. As long as the rotation remains commensurate and has 
a single wavevector, however, the induced density will have the same periodicity as that of the stripes. This can be 
seen as follows. The relation between the topological charge and the pseudospin texture is given by 



(r)) = ^-e abc S a (r) eijd % S b (r) d S c (r) . 



(50) 



where £y and e a bc are antisymmetric tensors and S (r) is a classical field with unit modulus representing the pseu- 
dospins. If we write a general solution as 
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S x (r) = sin [9 (r)] cos [<p (r)] , (51) 
S y (v) = sm[0(r)]sm^(r)], (52) 
S z (v) = cos[0(r)], (53) 

then the induced density takes the simple form 

Sn (r) = — sin [6 (r)] [V(^ (r) x V0 (r)] • z. (54) 

47T 

Writing a general solution for <p (r) in the form 

<p(r) = £(r) + Q-r, (55) 

the induced density can be written as 

Sn (r) = — sin [6 fr)l [VS (r) x V6» (r) + Q x V6> (r)l • z. (56) 

4-7T 



In the CSP and ISP, tp (r) = and so 



<5n (r) = — sin [9 (r)l [Q x V0 (r)l • 2. (57) 

47T 



There is no induced charge for stripes aligned with the parallel field [i.e. Q||V0(r)]. For other orientations of the 
stripes, the induce charge has the same periodicity as that of the angle (r) that defines the staggered magnetization 
in the xz plane. 

In our Hartree-Fock calculation, we find that the most favorable configuration is that of thp-stripes aligned with the 
parallel magnetic field. The same conclusion was also reported previously by Demler et alJi-3. Thus, we will choose 
a = in the rest of this paper. 



D. Incommensurate Stripe Phase (ISP) 

In analogy with the IL state, we define an incommensurate stripe phase (ISP) by the absence of rotation of the spins 
in the xy plane. This state is most naturally described by the original order parameters (p,y(q)) which, together 
with Gl,r{q, w n ), are non-vanishing for q € {G}- In particular, Gl^r (q — Q, w n ) = for any q € {G} so that 
the tunneling term in Eqs. (^)-(jl^) is zero. The order parameters of the ISP are thus found by solving directly 
Eqs. @-@ with t(q) = (or, alternatively, Eqs. @-(@) with t N (Q) = and Q = 0). While this is not an 
exact solution to the Hartree-Fock equations, what is ignored are the remnants of the solitons of the SLS which are 
relatively densejand leave only a very weak modulation of the pseudospin, except quite close to the transition out of 
the CSP stateEj. Setting t (q) = effectively ignores this weak modulation, and introduces a great simplification in 
the numerics. The resulting energy one finds for the ISP does not depend on the parallel magnetic field, just as for 
the IL state. We note that, because of the absence of tunneling in the Hartree-Fock equations, our approximate form 
of the ISP is equivalent to the coherent stripe state studied in Refs. and 0. 



V. PHASE DIAGRAM IN LANDAU LEVEL N = 2 



In the Hartree-Fock approximation, the transition between the CL to the IL states occurs when Ecl = Eil or, 
equivalently, when 



1 

4 


X (0) - X (Q)~ 


e Q 2 e 2 ± /4 




LI 







to 



(e 2 /n£)' 



(58) 



In Fig. 1, we show the contour lines in the d—B^ plane separating the CL from the IL for several values of the tunneling 
parameter. (We assume a typical density of n = 1 x 10 11 cm -2 in our calculation and in all the graphs presented 
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in this paper). We remark that the oscillations of (Q) lead to a complex sequence of transitions between the two 
liquid states with two reentrant CL phases for N = 2. We do not analyse these transitions in further detail, however, 
because at a value of d > d c (B.), the liquid states are higher in energy than the stripe states as indicated in Fig. 1, 
and the value of d c (B^ ) is lower than the value of d at which the reentrant CL states appear. It is nevertheless an 
interesting possibility that these reentrant transitions might occur when Landau level mixing or finite well thickness 
effects are included in the calculation. 

When stripe phases are taken into account, the phase diagram becomes even more complex because the oscillations 
of 1 7v (Q) also produce reentrant CSP and ISP states. Figure 5 shows the variation of the energy of the four states 
with respect to d for some typical cases. The interchange of the ISP and CSP is clearly visible in Fig. 5(a). The line 
called d c (i^ ) in Fig. 1 indicates the transition from one of the liquid states to one of the stripe states (either CSP 
or ISP). At small value of to, d c is almost independent of B . At larger values of d or B where tjy (Q) — > 0, the 
amplitude of the pseudo-magnetic field goes to zero and there is no gain in tunneling energy for the pseudospins to 
rotate. On the contrary, there is only an exchange energy cost. It is thus clear that the ground state must be the ISP 
at large d or B and the CSP only exists in a limited parameter region. Note however that, at large value of d, the 
LCR's are very narrow and the cost of the commensurate rotation in exchange energy is very small. The difference 
between the CSP and ISP energies in this limit is close to our numerical accuracy as can be seen in for large d in 
Fig. 5(b) and large in Fig. 6(a). 

From an experimental point of view, it is more interesting to study the behavior of the energy of the four states 
with respect to B at a fixed value of d. Referring to Fig. 1(a), we choose, for to = 0.01, d/£± = 1.4, a value greater 
than d c (B ) so that the ground state is one of the stripe states. We plot, in Fig. 6: (a) the energy of the CSP 
and ISP, (b) the variation in t^ {Q) , and (c) the corresponding variation of the order parameters (j) R>L (0)} C5P an d 
(Pr,l (°))jsp ■ ( We remark that, in the CL, (p RL (0)) = 0.5 and (p RL (0)) = while in the IL, (p RL (0)) = and 
(Prl (0)) = 0-5)- In the specific case studied in Fig. 6, the ground state goes through the sequence CSP-ISP-CSP- 
ISP-CSP-ISP when B is increased. A change in the sign of tjv (Q) leads to a cusp in the curve of the energy of the 
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CSP and to a corresponding change in sign in (p R L (0)). With our choice of phase, the pseudospins are oriented along 
x in the rotating frame for the CSP. When (Q) changes sign, the corresponding effective magnetic field is reversed 
and the pseudospins then point along —x. In computing Fig. 6, we have included only the CSP and ISP. But, we 
know from the study of the effect of B,, on the ground state at v = 1 that the commensurate state does not go directly 
into the incommensurate state. At a certain B,, = B..cl—SLS) which is smaller than the value B^ = B^^cl-il at 
which the CL — IL transition occurs, the commensurate state becomes unstable towards the formation of solitons. 
For B|| > B» t cL-SLSi the ground state contains a finite density of these defects. They form a soliton lattice whose 
lattice spacing depends on B.. . In analogy with the case v = 1, we expect that the same scenario will repeat itself here 
whenever the CSP approaches the ISP from below or_above. A precursor of this instability is seen in the softening of 
the low-energy modes of the CSP as described belowo. 

There are many possible transitions between the four states kept in our analysis. In Fig. 1, we indicate all transitions 
between the CL and IL states, but not all transitions into commensurate or incommensurate stripe phases that occur 
for d > d c . In general, we expect transitions between the CSP and the ISP whenever the tunneling matrix clement 
becomes small, as occurs in the liquid states. Thus, just as we see d c occurs near the first CL-IL transition in Fig. 1, 
we expect higher order CSP-ISP transitions near what are shown as the reentrant CL-IL transitions. Clearly, this 
system supports a remarkably rich phase diagram. 

How might one detect these phase transitions in experiment? One possibility, as we discuss below, is to study the 
behavior of the tunneling current with the parallel magnetic field. As we shall see this affords a clear signature of 
the transition between liquid and stripe phases. To compute the tunneling current, we need to know the collective 
excitations and response functions related to the density and pseudospin orders. We address this question in the next 
section. 



VI. RESPONSE FUNCTIONS IN THE TIME-DEPENDENT HARTREE-FOCK APPROXIMATION 



As mentioned above, the collective modes are important in determining the tunneling current in the double layer 
system in response to an interlayer bias potential. Their dispersions can also be measured experimentally via inelastic 
light scattering. To derive them, we compute the retarded density and pseudospin response functions, and track their 
poles. The response functions are obtained by analytical continuation of the two-particle Matsubara Green's functions 



Xi,j,k,i (q, q'; T ) = ~ N <t> ( T Pij (q, r ) Pk,i (-q'> °)) + N * (Pij (q)) (Pk,i (-q')> ■ 



(59) 



The response functions Xi,jkl (q^q'; 7 ") are computed in the time-dependent Hartree-Fock approximation. The cal- 
culation we do here is a generalization of that presented in Refs. || and [l3|. To include the parallel field and so the 
commensurate rotation, it is necessary to derive the equation of motion for the matrix of response functions defined 
in the rotating frame by 



P(q,qVO n ) 



( i*x R ,R,R t R (q, q') 7 i*xr, r ,l,r (q, q' - Q) i*Xr,r,r,l (q> q' + Q) i*x r , r .l,l (q,q')7* \ 

x r ,l,r,r (q - Q.q') 7 xr,l,l,r (q - Q>q'- Q) x r ,l,r,l (q - Q, q' + Q) x r , l ,l,l (q - Q>q') 7* 

xl,r,r,r{<i + QM)i xl,r,l,r(<i + Q,<i' - Q) x l ,r,r,l (q + Q, q' + Q) XL,i?,L,L(q + Q,q')7* 

V iXL,L t R,R (q> q') 7 7Xl,l,l,_r (q, q' - Q) 7Xl.l„r.l (q> q' + Q) ixl.l.l.l (q,q')7* / 

(60) 



where 7 on the left(right) of P(q, q';Q n ) stands for 7(q)(7(q')) an d is a Matsubara bosonic frequency. In the 
time-dependent Hartree-Fock approximation (TDHFA) , the summation of ladders and bubble diagrams gives 



P (q, q'; ifi„) = P° (q, q'; iO n ) + \ £ P° (q, q"; ifi„) Y (q") P (q", q'; 



(61) 



q" 



The vertex corrections involve the matrix of interactions 



e 

kI± 



( H (q) - X (q) 



V i?(q)[7(q)] 2 



^0 
-X- (q) 





H( q )[r(q)f\ 

~° 
-X+ (q) 

H(q)-X(q)J 



(62) 



where we have defined 



X±(q)=*(q±Q). 



(63) 
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FIG. 6: Study of the CSP-ISP transition at N = 2 for to = 0.01 and d/£± — 1.4 showing (a) the variation of the energy of these 
two phases with B\\ , (b) the variation of the tunneling amplitude tM (Q) with Bu and (c) the variations of the order parameters 

(Pr,l (0)) CS p and (Pr,l (°))i SP with B \\- 



Using the Hartree-Fock Hamiltonian, Hhf, written in terms of the (/jj j (q)), we obtain an equation of motion for 

the Hartree-Fock response functions P° (q, q'; iVt n ) using -§- (...) = [H^f — I^N, (...)]. Combining with Eq. (|l|), we 
obtain after a lengthy calculation an equation of motion for the TDHFA response functions than can be put in the 
simple matrix form 

53 ^ [iMlJ^S^,, - F i>k (q, q")] P k ,j (q", q') = Aj (q, q') ■ (64) 

k q" 

For completeness, we give the definitions of the matrices F and D in the Appendix. Notice that the dynamical 
response can be computed from a knowledge of the static order parameters {(Pi j (q))} alone. 
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A. Response Functions of the Commensurate Liquid 



For the CL, we can solve analytically Eq. (|64|) to get (after the analytical continuation iQ„ — ► u + iS and with 
( = sign[t N (Q)]) 

/ xr, r ,r,r (q, q') Xfl,i?,L,_R(q 5 q' - Q) xr, r ,r,l (q> q' + Q) x r , r ,l,l (q, q') \ 
xr,l,r, r (q - Q, q') Xii,L,i,fi(q-Q>q'-Q) x;? : L„R,L(q-Q,q' + Q) xr,l,l,l (q- Q>q') 
xl, r ,r,r (q + Q, q') Xi,fl,i,fl(q + Q,q'-Q) x l ,r,r,l (q + Q> q' + Q) x l ,r,l,l (q + Q,q') 
V xl,l,r, r (q. q') Xi,i,i,fl(q 5 q'-Q) xl,l,.r,l (q, q' + Q) x l ,l,l,l (q, q') / 

/ _(6 + c) -ft(u + i<5) 7 fi(w + i(5)7 -0 + c) 7 2 ^ 



(65) 



(C/2)<5c 



(a; + i<5) 2 — w 2 



CL 



(q) 



— H (uj + i5) 7* 
H{uj + iS) 7* 



2 Re (a) 
-2 Re (a) 



\-(7*) (6 + c) H(lu + i5)Y 



-2 Re (a) 
2 Re (a) 
-H (u> + «(5) 7* 



h (u> + iS) 7 
— H (ui + iS) 7 
+ (b + c) J 



where 



a(q) 
c(q) 



2|*jv (Q)| 
(Q)\ 



e 

2 



X (Q) - (g) 



2|MQ)I+ ~H -*-(«) 



(66) 
(67) 
(68) 



There is a unique collective mode corresponding to an elliptical motion of the pseudospins about the x axis, in the 
yz plane, with a frequency given by 



Hlu cl (q) = y/[a (q) + a* (q)] [b (q) + c (q) 



(69) 



(We remark that, for q = 0, Jf (q) — ZL(q) cos [q x Q/jJ — » —d/lj_). Note w^l (q) tends to a finite value for any 
tjv(Q) > in the long- wavelength limitlla. 



B. Response Functions of the Incommensurate Liquid 



For the IL, the response functions are obtained by setting tjy (Q) — and Q — 0. We get 

/ Xfl,fl,fl,ji(q,q') xr, r ,l,r (q. q') x^ii^Cq.q') x r , r ,l,l (q, q') \ 

XH,i,ie 1 ii(q 1 q') xr,l,l,r (q, q') XR,i,K,i(q,q') xr,l,l,l (q, q r 

xz,,ij,H,ij(q»q') x l ,r,l,r (q> q') xl, r ,r,l (q> q') x l ,r,l,l (q>q' 

V x l ,l,r,r (q> q') x l ,l,l,r (q> q') x l ,l, r ,l (q> q') Xi.i.i.iC^q') 7 

25 -ft(w + i<5) fi(w + i<5) 
(1/2) <5 q ,q/ ( -ft( w + i<5) 24 -24 



where 



(w + itf) -^| L (q) 



A(q) = - 



(w + i<5) 
-25 



-24 
h(uj + i(5) 



24 



e 



5(q) = = -7- 



e 



ff(q) 

X(0) 



■if(q)+X(0)-X(q) 
*(q)" 



The collective mode is gapless and given by 

fiw JL (q) = 2^4 (q)5(q). 
To linear order in g (and for ^ 0) , this becomes 



fiw /L (q) 



1\ <9 2 X(fc) 



<% 2 



X(0)-X(Q) 



Thus for long wavelengths the dispersion is linear in q and gapless even if t > 0. 



-25 
H(lu + iS) 
-H (to + i<5) 
25 



(70) 



(71) 
(72) 

(73) 
(74) 
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C. Response Functions in the Commensurate Stripe Phase 

We study the dispersion relation of the lower-energy collective modes of the CSP by a numerical solution of the 
systertLof equations given by Eq. ([m]). In the absence of tunneling and parallel magnetic field, there are two Goldstone 
modesu. One is a phonon mode with an in-phase motion of the density of the stripes in the two wells accompanied 
by an x — z motion of the pseudospins. The second is a pseudospin wave mode with an out-of-phase motion of the 
densities and a small rotation of the pseudospins about their equilibrium position. As explained in Ref. ^, the low- 
energy modes can be described by a pseudospin-wave theory where the effective periodicity is set by the separation 
between neighboring linear coherent regions (LCR's) which is half the separation between the stripes in a given layer. 
If we take the the Brillouin zone edges (in the x direction) to be ±2n/a instead of ±7r/a, then the phonon and 
pseudospin wave mode are both part of the same low-energy mode, the former dispersing from the zone center and 
the latter from the zone edge. This is very clearly seen from Fig. 2(a) and Fig. 2(b) where we show the full dispersion 
relation of these modes in the ISP and in the CSP with Bn = OT. In Fig. 2(b), we see that the gap at the zone 
edge (and so in the pseudospin wave mode) is strongly increased from its bare value 2to by self-energy and vertex 
corrections. The phonon, by contrast, has vanishing energy at the zone center. For Bu = 0, the phonon mode has 
lu (fc||, k± = 0) ~ k 2 and lu (fey = 0, fc_i_) ~ k± and is independent of to. For B\\ — and to = 0, the pseudospin wave 

mode has lu (fc||, fc^ = 0) ~ k\\, lu (k\\ = 0, fc_i_) ~ k± (see Ref. [)]). 

In Fig. 6, we showed the behavior of the renormalized tunneling term, t/v (Q), as B is increased at d/£± = 1.4 for 
N = 2. We now look at the corresponding behavior of the dispersion relation for the phonon and pseudospin wave 
modes. (As discussed above, these two modes are both part of the same collective mode in an extended Brillouin 
zone such as that used in Fig. 7). Figure 7 shows the dispersion relations of these collective excitations as B» is 
increased until slightly above the first value where tN (Q) — where the system is in the ISP. The parallel field has 
an important effect on the lowest-energy modes. The gap in the pseudospin mode (the part of the dispersion near the 
zone edge in Fig. 7) is progressively closed as t N (Q) — > 0. But, these modes do not become completely dispersionless 
along fen = 0. Before the CSP-ISP transition, these low-energy modes become unstable at the zone center at some 
value B* < B^ CSP _ ISP . The instability persists in a small region around B^ CSP _ ISP . (For the case represented 
in Fig. 6, we found that the first instability region is from Bu = 0.34 to Bu = 0.40). Above B^ CSP _ ISP , the ground 
state is stable again and the gap increases and decreases again until the next zero of tjv (Q) = where the same 
scenario is repeated. Finally, just before the final transition from the CSP to the ISP, i.e. at B.. =1.6 for the case 
represented in Fig. 6, the CSP becomes unstable again. 

In analogy with the case of filling factor v = 1, we expect that these instabilities are precursors of a transition from 
the CSP to a soliton lattice state. For v = 1, the instability occurs for B* > B^ CL _ SL g and the transition from 
the commensurate liquid to the SLS is first order. We do not know if this is the case for the CSP to ISP transition 
as we could not calculate B^csp-sls with our Hartree-Fock method. In any case, these instabilities need not be 
worried about when considering the low energy spectrum: the introduction of solitons into the groundstate will open 
gaps in the collective mode spectrum at various values of but should not provide any new gapless modes since 
translational symmetry has already been broken by the stripes. 

Figure 8 shows the dispersion relation of the collective modes in the direction of the stripes for several values of the 
parallel magnetic field. It is clear that a change in B.. has a profound impact on the long-wavelength dispersion of the 
phase mode (lowest energy gapped mode at k^ = 0) as well as on the roton minimum. In constrast, the higher-energy 
modes are much less affected by the parallel field. (This is also true for dispersion perpendicular to the stripes) . The 
phase mode becomes nearly linear in fcii at k± = near the zeros of tjv (Q)- The long- wavelength phonon dispersion 
changes slowly with B^ . At small B^ , we find lu (k\\ — > 0, k± = 0) ~ k^ while at B^ = 0.8 T (and for d/£± = 1.4), the 
dispersion crosses over from quadratic to linear in fc|| at a relatively small value of fc|| [see Fig. 8(d)]. 

As the separation between the wells increases, both the CSP and ISP are unstable and the true ground state is that 
of a modulated stripe state described in Ref. ^. For to/ (e 2 /n£) = 0.01, we find that the critical spacing is d/£± ps 2.0. 

VII. TUNNELING RESPONSE IN THE COMMENSURATE STRIPE STATE 

In this section we compute the interlaycr tunneling I — V to second order in the tunneling matrix element io- 
Since to is assumed to be small, the commensurate states can only exist in a very narrow range of parameters d and 
B.. , so that they may effectively be ignored. We are thus interested in the IL-ISP_txansition, and whether there are 
distinguishable signatures of these two phases in the tunneling I — V . Recent workal!3 has demonstrated that a peak 
in the tunneling I — V as a function of B.. maps out the Goldstone mode associated with spontaneous interlayer 
coherence in the IL. As we now show, this result generalizes to the ISP, and the behavior of the peak with B.. clearly 
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FIG. 7: Dispersion relation to (k\\ —0,k±) of the lowest-energy collective mode in the CSP for N = 2, to/ (e 2 / ' K,t) — 0.01, d/£± = 
1.4 and several values of the parallel magnetic field. The [2n/a, An /a] region can be translated left by 2n/a to give the pseudospin 
wave dispersion. 



distinguishes the two phases. I— . 

Our formalism closely follows that of MaharO. The tunneling current is defined by 



J(t) = -e- 



dNL 

dt 



dN R 

1 dt ' 



(75) 



where Nl^^ = J^x c x L(R) Cx , L (R) = ^<f>PL(R) L(R) = 0) i s the number of electrons in the left (right) well. I (t) 
represents the interlayer current in the presence of a bias voltage while the tunneling operator connecting the two 
layers is adiabatically switched on, reaching its full value at time t — 0. The complete Hamiltonian of the system, 
including a bias between the wells, is written as 



H = H' + H T = H a - [i L N L - fi R N R + H T , 



(76) 



where eV = \l l — [i R is the potential bias and Ht is the tunneling hamiltonian given in Eq. (??). Written in terms 
of the p operators, it takes the form (Q —Qx) 



H T = -N^t N (Q) [p R>L (-Q) + p L<R (Q)] 



(77) 



The term Hq in H contains the Coulomb interaction and kinetic energies. Since Ht is the only term that can change 
particle numbers in one well, we have 



(/.V/,. 
dt 



= i [H T , N r ] = -N^t N (Q) [p R>L (-Q) - p L ^ R (Q)] 



(78) 



In all four states considered in this paper. (p R L (— Q)) = (pl R (Q)) so that there is no current in the ground state. 
Using standard linear response theoryo, the current is given by 



(Ht)) 



dt' 



H T (t'),I(t) 



H' 



(79) 
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FIG. 8: Evolution of the dispersion relations uo (k\\ , k± = OJ of the lowest-energy collective modes of the CSP as a function of 
the parallel magnetic field for N = 2, to/ (e 2 / ' kI) = 0.01 and d/tj_ = 1.4 

where the time dependence of the operators is given by I (t) = e iH '*/ fi 7 (0) e - lff '*/ fi . After some algebra, we find the 
familiar result 

I (t) = ~2^iV \t N (Q)| 2 Im [ Xr , l<l ,r (-Q, -Q,fiw = eV)] (80) 
-2^iV |i w (Q)| Im [e- 2ieW Xjt , L , fl , L (-Q, +Q,fio;=eF)] . 

The response functions in Eq. ( p0| ) must be evaluated in a ground state where tunneling and so the parallel magnetic 
field are absent. Thus, only the response functions of the IL or ISP are relevant in this approach to the tunneling 
I — V. In the liquid state where we have the analytical result given by Eq. (|70|), the second term in Eq. (|o|) does not 
contribute to the tunneling current and we find we the tunneling current of the IL phase for small to, 



Iil (t) = ^ \t N (Q)\ 2 *J §||<5 (eV - fko IL (Q)) , (81) 



where A (Q) and B [Q) are given by Eqs. ( J7l| ) and ( [72D above and lojl (Q) is the Goldstone mode of the IL state 
whose dispersion is given by Eq. ( |73| ) . This result is identical to that of Ref . [h], and as pointed out there the peak 
in the tunneling current can be used to follow the dispersion relation of the Goldstone mode. In the experiment of 
Spielman et al.B, the tunneling amplitude is extremely small and the expression above can be applied provided one 
introduces a phcnomcnological broadening due to disorder and thermal fluctuations. This can be achieved rather 



naturally by substituting a Lorentzian for the delta function in Eq. (81), although the microscopic mechanism leading 
to this broadening is currently not understood. 
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Whereas in the IL state, the tunneling current peaks near zero bias only for B.. = 0, in the ISP the Goldstone 
mode folds back upon itself, continuously deforming into the phonon mode as the Brillouin zone edge is approached. 
This means the peak in the tunneling will return to zero bias when B.. reaches a value such that Q = ±2mr/a, with 
n = 0, 1, 2, 3, . . . provided ijv (Q) is not zero. This is illustrated in Fig. 3, which is produced using the first (normal) 
term in Eq. ( |80| ) with a phenomenological broadening parameter 5. This non-monotonic behavior of the tunneling 
current peak is an unambiguous signature of the ISP that distinquishes it from the IL. In this way one may probe 
the transition from IL to ISP via mcasiu'erncnt of the interlaycr tunneling current. Note that the tunneling current 
vanishes at large B.. where tjy (Q) — > OEa. 

The perturbative approach used in this section to compute the tunneling current does not capture the full complexity 
of the phase diagram discussed in Sec. V. To probe the various phase transitions discussed above and the corresponding 
change in the collective excitations, the response functions must be evaluated in a ground state that includes the 
parallel magnetic field. That is, the response functions must be evaluated non-perturbatively with respect to the 
tunneling matrix element i - Work on this problem is currently in progress. 



We have shown that, in higher Landau levels (N > 1), a parallel magnetic field can induce a complex sequence of 
phase transitions in the two-dimensional electron gas of a double-quantum-well system. When the separation between 
the wells or the parallel magnetic field are varied, transitions occur between commensurate and incommensurate 
liquid or stripe states. Working in the generalized random-phase approximation, we have computed the spectrum 
of collective excitations of these four ground states showing how the pseudospin-wave mode of the incommensurate 
liquid phase is modified by the formation of the stripes in each well and by the parallel magnetic field. We believe 
that the (incommensurate) stripe state could be distinguished from the (incommensurate) liquid state in a tunneling 
experiment because of the distinctive signature due to the periodicity of the ISP. 
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APPENDIX A: DEFINITION OF THE MATRICES F AND D 



With 



S. 



< qj q' = sm 



q x q'l 
2~ 




(Al) 



The matrices F and D of Eq. (p4) are given by 




(A2) 




(A3) 
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A,3 (q,q') = 



P^q-q')^ 




P^ (q - q')) 7 q , q < - (P H ' H (q - q'); 7 q , q 

-(p R ' L (q-q'))7q, q ' 



(A4) 



\ 



p fl ' L (q-q'))7 q , q ^ 



p L ' fl (q-q')>7 q , q < 



A.4 (q, q') = 

^ 2iSq M > (p^ L (q- q')) ^ 
and (in these definitions, the tunneling {Q) and all matrix elements are in units of e 2 /n£) 

Fi.i (q, q') = -2*S q , q , [[H (q - q') - X (q - q') - if (q') + X (q')] (p R ' R (q - q') 
-2iS q , q ,£ (q - q') (7* (q - q')) 2 (p^ (q - q')) - 



(A5) 



(A6) 



f2,i(q,q') - -t N (Q)S^- H(q')-A:(q')+X-(q-q') (p fl ' L (q - q')) 7q , 



<q,q 



ff(q')7^(q') (p fl ' L (q-q'))7 p , p ' 



(A7) 



F 3 ,i (q,q') = tN (Q) <5 q , q ' + \ H (q') — X (q') + X + (q - q') (p (q - q') ) 7 q , q < 
^(q')7 2 (q')l (p L>fl (q-q')W q <, 



F 4 ,i(q ) q , ) = 2iS q , q , H (q') 7 2 (q') p^(q-q')), 



i r i,2(q,q , ) = -^(Q)*q,q'- X+(q-q')-X_(q') ( p (q — q') / 7q,q' i 



^2,2 (q,q') = - 



if (q - q') - X (q - q') + X_ (q')J (q - q')) 7 q ,c 

+i? (q - q') ( 7 (q - q'))* (p*'* (q - q')) 7 q , q < 



if (q - q') - X (q - q') + X_ (q')j (p^ (q - q')) 7 q , q < 
H (q - q') ( 7 * (q - q')) 2 (p L < L (q - q')) 7 q , q <, 



(A8) 
(A9) 
(A10) 
(All) 



F 3 .2 (q,q') = 0, 



^4,2 (q, q') = t N (Q) * q , q , + (J^J [X+ (q - q') - X_ (q')] (p L ^ (q - q')) 7q , q „ 



(A12) 
(A13) 



J , i,3(q,q')=*iv(0)*q,q'+ (^-) X-(q-q')-^+(q') (p ' (q " q') ) 7 q , q ', 



^2,3 (q,q') = 0, 



(A14) 
(A15) 



*3,3(q,q') = - ff(q-qO-X(q-q')+X + (qO ( ~p L ' L (q - q') ) 7* 



'q,q 



+H (q - q') (7* (q - q')) 2 (~p L L (q - q')) 7 q , q ' 



H (q - q') - X (q - q') + X + (q')j (q - q')) 7 q , q ' 

H (q - q') (7 (q - q'))' (p*'* (q " q') ) 7o, q < , 
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(A16) 



^4,3 (q, q') = -t N (Q) <W - X- (q - q') - X+ (q') (p K ^ (q - q') ) 7q , 



i>q'> 



F 1A (q, q') = 2zS q , q , (q - q') ) H (q') (7* (q')) , 



F 2 , 4 (q,q') = MQ)<W + jL(q-q') + ff(q')-*(q') (?*' £ (q - q') ) 7 q , q < 



ff(q')(7* (q')) 2 



p (q-qj/7 qq ' 



^3,4 (q, q') = -t N (Q) <W - X+ (q - q') + if (q') - X (q') (q - q') ) 7 q>q , 



H(qO(7*(q')r (p L '*(q-q')>7 q 



^4,4 (q,q') = -2i5 q , q ,ff(q- q ')7 2 (q-q')(p^(q-q'; / 

-2i£ q , q , [if (q - q') - X (q - q') - if (q') + X (q')] (p L ' L (q - q'; 



(A17) 
(A18) 

(A19) 

(A20) 
(A21) 
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